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ABSTRACT. Calving activity at the front of tidewater glaciers is characterized by a large 
variability in iceberg sizes and inter-event intervals. We present calving-event data obtained 
from continuous observations of the fronts of two tidewater glaciers on Svalbard, and show 
that the distributions of event sizes and inter-event intervals can be reproduced by a simple 
calving model focusing on the mutual interplay between calving and the destabilization of 
the glacier front. The event-size distributions of both the field and the model data extend 
over several orders of magnitude and resemble power laws. The distributions of inter-event 
intervals are broad, but have a less pronounced tail. In the model, the width of the size distri- 
bution increases with the calving susceptibility of the glacier front, a parameter measuring 
the effect of calving on the stress in the local neighborhood of the calving region. Inter- 
event interval distributions, in contrast, are insensitive to the calving susceptibility. Above 
a critical susceptibility, small perturbations of the glacier result in ongoing self-sustained 
calving activity. The model suggests that the shape of the event-size distribution of a glacier 
is informative about its proximity to this transition point. Observations of rapid glacier 
retreats can be explained by supercritical self-sustained calving. 
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Iceberg calving plays a key role in glacier dynamics and, 
hence, in how tidewater glaciers and ice sheets respond to 
climate change, thereby impacting predictions of sea level 
rise in the future (van der Veen, 1997; O'Neel et al., 2003; 
Benn et al., 2007a; Nick et al., 2009). So far, the mechanisms 
underlying the calving dynamics are only partly understood. 
To summarize the potential controls affecting iceberg calving, 
Benn et al. (2007a) proposed the following classification: first- 
order controls determining the position of the glacier front, 
second-order controls responsible for the calving of individ- 
ual icebergs, and third-order controls related to the calving 
of submarine icebergs. The first-order control on calving is 
the strain rate resulting from spatial variations in the glacier 
velocity, responsible for the opening of crevasses. Crevasse 
formation is reinforced by the presence of liquid water, either 
from surface melt or rain events. Second-order controls are 
processes weakening the glacier front and favoring fractures, 
like the presence of force imbalances at the glacier front re- 
sulting from the margin geometry, undercutting of ice and 
torque due to buoyancy. Third-order controls trigger subma- 
rine iceberg calving by processes like the formation of basal 
crevasses, tides and buoyancy. 

The majority of previous studies describe calving by means 
of macroscopic measures like the overall calving rate (calv- 
ing speed), i.e. the total ice loss at the glacier front within 
rather long time intervals. Corresponding models ("calving 
laws") relate the dynamics of the calving speed to parame- 
ters like the water depth (Brown et al., 1982; Oerlemans et 
al., 2011), the height-above-buoyancy (van der Veen, 1996; 
Vieli et al., 2001, 2002), the penetration of surface and basal 



crevasses arising from the longitudinal strain near the calving 
front and surface melt (Benn et al., 2007b; Nick et al., 2010; 
Otero et al., 2010), and more general glacier characteristics 
like the ice thickness, the thickness gradient, the strain rate, 
the mass-balance rate, and backward melting of the terminus 
(Amundson & Truffer, 2010). 

So far, only few studies have been dedicated to a description 
of calving dynamics at the level of individual calving events. 
Continuous monitoring of individual events directly at the 
glacier terminus is difficult. In consequence, data are sparse. 
Previous studies of the calving-event statistics are based on 
occasional or non-continuous observations of calving events 
(Washburn, 1936; Warren et al., 1995; O'Neel et al., 2003), 
or on indirect measurements, e.g. event sizes obtained from 
icebergs floating in the sea (Budd et al., 1980; Orheim , 1985; 
Wadhams , 1988) or from seismic activity (O'Neel et al., 
2010). The available data indicate that event sizes are highly 
variable and broadly distributed (e.g. Bahr, 1995; O'Neel et 
al., 2010). However, distributions of event sizes obtained from 
floating icebergs in the sea are likely to be biased due to melt- 
ing and disintegration (Neshyba, 1980; Bahr, 1995). Seismic 
measurements can reliably detect only large calving events 
(Kohler et al., 2012). In addition, estimating calving-event 
sizes (ice volume) from seismic-event magnitudes is problem- 
atic unless the relation between these two measures is clearly 
established (e.g. through calibration by means of visual obser- 
vations; see Kohler et al., 2012). Hence, the resulting size dis- 
tributions may be distorted. In principle, calving activity at 
the single-event scale could be monitored by means of repeat 
photography (O'Neel et al., 2003), laser scanning, or ground- 
based radar (Chapuis et al., 2010). So far, however, no such 
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data has been published yet. Here, we present single-event 
data obtained from continuous visual observations directly 
at the fronts of two tidewater glaciers on Svalbard. Our data 
confirm that both the sizes of individual calving events and 
the time intervals between consecutive events are broadly dis- 
tributed. The mechanisms underlying this calving variability 
are so far unknown. It is unclear whether it reflects variabil- 
ity in external conditions, e.g. fluctuations in temperature 
or tides, or whether it is generated by the internal calving 
dynamics itself. 

As fluctuations in external conditions can hardly be con- 
trolled in nature, disentangling these two potential sources of 
variability requires a model of the calving process. A descrip- 
tion of size and timing of individual calving events is beyond 
reach of the macroscopic continuum models focusing on the 
overall calving rate (see above). Bahr (1995), Bassis (2010) 
and Amundson & Truffer (2010) proposed models accounting 
for the discreteness of calving. In the model of Amundson 
& Truffer (2010), calving events are triggered when the ter- 
minus thickness decreases to some critical value. According 
to this model, the event sizes and inter-event intervals are 
fixed for given model parameters. Variability in these mea- 
sures can result from fluctuations (e.g. seasonal variations) 
in these parameters. Bassis (2010) describes the motion of 
the glacier terminus in one dimension as a stochastic process. 
In the study of Bahr (1995), calving is modeled as a percola- 
tion process in a two-dimensional lattice representing a region 
close to the glacier terminus. In this model, microfractures 
are randomly and independently generated according to some 
cracking probability. Calving events occur whenever a section 
of ice is surrounded by a cluster of connected microfractures. 
In this model, calving in one region of the model glacier has 
no effect on the state of the rest of the glacier. Both the model 
of Bahr (1995) and the one of Bassis (2010) are stochastic. In 
our study, we propose an alternative calving model focusing 
on the mutual interplay between calving and the destabiliza- 
tion of the local neighborhood of the calving region. Although 
the calving dynamics of our model is fully deterministic, it 
generates broad distributions of event sizes and inter-event 
interval distributions which are consistent with the field data, 
even under stationary conditions. 

The article is organized as follows: After describing the 
acquisition of the field data, the calving model and the data 
analysis (Sec. 2), we present the event-size and interval dis- 
tributions obtained from the field data and show that they 
are reproduced by the calving model (Sec. 3.1). The model 
predicts a large variability in the calving process even under 
stationary external conditions. This is supported by the field 
data showing that the shape of the size and interval distri- 
butions is not affected by registered fluctuations in climatic 
conditions (Sec. 3.2). Further, we discuss a prediction of the 
model which may be of significance for diagnosing the sta- 
bility of a glacier: the model glacier exhibits a critical point 
at which it enters a regime of ongoing self-sustained calving 
(Sec. 3.3), a regime which may be related to observations of 
rapid glacier retreats. Finally, we demonstrate that the calv- 
ing model is consistent with the field data in the sense that the 
size of future calving events is barely predictable from past 
events (Sec. 3.4). In the last section of this article (Sec. 4), we 
discuss the conclusions and consequences of our work, embed 



the results into the literature and point out limitations and 
possible extensions. 

2. METHODS 

2.1. Acquisition of field data 

Calving activity was monitored at two tidewater glaciers on 
Svalbard: Kronebreen and Sveabreen. In this region, midnight 
sun lasts from April 18 to August 24, thereby allowing for a 
continuous (day and night) monitoring of calving activity. 
Kronebreen is a grounded, polythermal tidewater glacier, lo- 
cated at 78°53iV, 12°30-E, approximately 14 km south-east 
of Ny-Alesund, western Spitsbergen (Fig. 1A). Kronebreen is 
one of the fastest tidewater glaciers on Svalbard with an aver- 
age front velocity between 2.5 and 3.5 m/d during the summer 
months (Rolstad & Norland, 2009) . At the end of August 
2008, the terminal ice cliff had an elevation ranging from 5 
to 60 m above the fjord surface (Chapuis et al., 2010) . Four 
persons observed the calving front of Kronebreen during a 
total of 16 days, split into two periods, from August 26th, 
2008, 19:00, to September 1st, 2008, 05:11, and from August 
14th, 2009, 00:00, to August 26th, 2009, 16:00 (GMT). The 
observation camp was located approximately 1.5 km west of 
the glacier front (Fig. 1, open white triangle). About 90% 
of the front was visible for the observers. The second glacier, 
Sveabreen, is a 30 km long, grounded tidewater glacier located 
at 78°33A'', 14°20£, flowing in the northern part of Isfjorden 
(Fig. IB). The observation of Sveabreen was part of a Youth 
Expedition program. About 45 persons observed the front of 
Sveabreen during 4 days from July 17, 2010, 14:40, until July 
21, 2010, 15:00 (GMT). The camp was located approximately 
500-700 m from the glacier front and offered a good and open 
view of the front. 

For each calving event, we registered the type, location, 
time and size. For our analysis, we pooled all event types 
(avalanche, block slump, column drop, column rotation, sub- 
marine event; see O'Neel et al., 2003; Chapuis et al., 2010). 
The event times were monitored with a relative accuracy of 
10 seconds. Due to delays between the actual events and 
the hearing and registering, we assigned an absolute tem- 
poral precision of 1 minute to the inter-event intervals r 
(time between two consecutive events). Following the semi- 
quantitative approach introduced by Warren et al. (1995), we 
monitored the perceived size ip £ {1, 2, . . .} of each calving 
iceberg on an integer scale (O'Neel et al., 2003). The small- 
est observable events (ip = 1) correspond to icebergs with a 
volume of about 10m 3 , the largest (ip = 11) to more than 
10 5 m 3 (collapse of about l/5th of the glacier front). During 
common observation periods, the perceived event sizes ip reg- 
istered by different observers could be compared. Based on 
these data, the error (variability) in ip is estimated as 1. The 
perceived iceberg sizes ip were mapped to the actual iceberg 
volumes fi by means of photogrammetry: Repeat photographs 
were automatically taken in 3-second intervals from a fixed 
location (star in Fig. 1) using Harbotronics time-lapse cam- 
eras (see e.g. Chapuis et al., 2010). In the resulting photo 
material, we identified 18 calving events which were simulta- 
neously observed visually. The approximate volume of each 
event was obtained from the estimated iceberg dimensions, 
and compared to the perceived size ip (see Fig. 2). The rela- 
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Fig. 1: Aerial pictures of Kronebreen taken in August 2009 (A) and Sveabreen taken in August 2010 (B). Locations of camp 
and time-lapse camera are marked by a triangle and a star, respectively. Inset: map of Svalbard showing the location of the 
two glaciers. 



tion between /j, and ip is well fit by a power law (dashed line in 
Fig. 2; correlation coefficient c = 0.68 in double-logarithmic 
representation) : 

At = 12.6.^ 3 ' 87 . (I) 

Note that the empirical power-law model (1) is consistent 
with psychophysical findings ( "Stevens' power law" ; Stevens, 
1957). Using (1), we converted the perceived iceberg sizes ip 
for all visually monitored events to an estimated volume \x 
(in units of m 3 ). 

2.2. Calving model 

Breaking of ice, formation of fractures (crevasses) and, in 
turn, calving are a consequence of internal ice stress (Benn et 
al., 2007a). Several mechanisms contribute to the build-up of 
stress at the glacier front, e.g. glacier- velocity gradients, buoy- 
ancy, tides, or changes in the glacier-front geometry due to 
calving itself. In this study, we focus on the mutual interplay 
between calving and the destabilization of the front, mediated 
mainly by a loss of buttress in the neighborhood of the calv- 
ing region, but also by positive feedback loops between calv- 
ing, reduction in ice burden pressure, increase in buoyancy, 
terminus acceleration and, in turn, increase in calving rates 
(see e.g. Benn et al., 2007a, and references therein). We show 
that this interaction loop between calving and front desta- 
bilization alone is sufficient to explain the large variability 
of iceberg sizes and inter-event intervals observed in the field 
data. In our model, all other stress contributors are treated as 
"external stress" or described by parameters. Keeping these 
parameters constant allows us to study the glacier dynamics 
under (ideal) stationary external conditions. 

Overview. In our model, the glacier front is described as 
a two-dimensional, discretized rectangular plane, subdivided 
into cells (Fig. 3C). Each cell corresponds to a unit volume 
of ice. The state of a cell is characterized by its internal ice 
stress. If this stress exceeds a critical level, the ice breaks, the 



cell "calves" and its stress is reset to zero (Fig. 3 A). Calving 
of a cell increases the stress in neighboring cells (Fig. 3C), 
mainly as a consequence of a loss of buttress. Hence, ini- 
tial calving of individual cells can trigger calving avalanches 
involving larger regions of the glacier front. We probe the 
model glacier by applying small perturbations (small stress 
increments) to randomly selected, individual cells. The total 
number of cells participating in an avalanche triggered by 
a single perturbation defines the event size fi. The distance 
between two consecutive events of non-zero size corresponds 
to the inter-event interval r. In the following, the different 
model ingredients are described in detail. 

Model geometry. Our calving model focuses on the calv- 
ing dynamics at the glacier front (terminus). For simplicity, 
the front is described as a two-dimensional rectangular plane 
of width W and height H. The front is discretized, i.e. subdi- 
vided into W H cells with coordinates {x, y\x = 1, . . . , W; y = 
1, . . . , H} (y — and y = H correspond to the sea level and 
the height of the glacier terminus above sea level, respectively; 
see Fig. 3C). Each cell represents a unit volume of ice. Note 
that the model neglects the third spatial dimension perpen- 
dicular to the front plane. 

Stress dynamics and calving. The internal ice stress 
in a cell at position {xy} at time t is described by a scalar 
variable z X y(t) (Fig. 3A). The cell calves at time t l xy if its 
internal stress exceeds a critical value of .Zcrit — 1 ("yield 
stress"; see e.g. Benn et al., 2007a, and references therein), 
i.e. if z X y(t l kl ) > 1 (triangle in Fig. 3A). The cell's calving 
activity can be described mathematically as a sequence of 
calving times {t xy \i = 1,2,...}, or, more conveniently, as 
a sum of delta pulses, s X y(t) = £\ S(t — t xy ) (triangles in 
Fig. 3A,B). After the cell at position {xy} has calved, it is 
assumed to be replaced by a "fresh" cell representing ice in 
a deeper layer. In the model, this replacement is captured by 
instantaneously resetting the stress at position {xy} to zero 
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Fig. 2: Measured iceberg volume fi versus perceived iceberg size ip (log- log scale) for 18 individual calving events (symbols). 
Error bars depict estimated volume measurement errors. The dotted line represents the best power-law fit (1) (linear fit in 
log-log representation). 



(triangle in Fig. 3A). Note that the geometry of the model 
(see above) is not altered by calving. We assume that the 
dynamics of the internal ice stress z X y{t) represents a jump 
process which is driven by calving of neighboring cells and 
external perturbations (triangles in Fig. 3B). Mathematically, 
the (subthreshold, for z xy < 1) stress dynamics can be de- 
scribed by 



dZxyjt) 

dt 



W H 



= E E J kl S kl(t) + JcxtsZ(t) - S X y{t) . (2) 



k=l 1=1 



Here, the left-hand side denotes the change in stress at time 
t (temporal derivative). The right-hand side (rhs) of (2) de- 
scribes different types of inputs to the target cell {xy}. In 
the absence of these inputs (i.e. if the rhs is zero), the stress 
level z X y remains constant. The first term on the rhs corre- 
sponds to the stress build-up due to calving in neighboring 
cells: calving of cell {kl} at time t leads to an instantaneous 
jump in z X y with amplitude Jf, (Fig. 3A,B,C). The second 
term represents stress increments as a result of external per- 
turbations s^ t (t) with amplitude J ox t- For simplicity, we as- 
sume that these external perturbations are punctual events 
in time (delta pulses), i.e. s^ t (t) = J2i^(t — t l cxt xy ). The 
last term on the rhs of (2) captures the stress reset after 
calving of cell {xy} (as described above) and is treated as a 
negative input here. Note that the single-cell calving model 
described here is identical to the "perfect integrate-and-fire 
model" which is used in a variety of other fields studying 
systems of pulse-coupled threshold elements like, for example, 
networks of nerve cells (Lapicque, 1907; Tuckwell, 1988), sand 
piles (Bak et al., 1987, 1988), or to investigate the dynamics 
of earthquakes (Herz et al., 1995). 

Interactions between cells. Calving of a cell at posi- 
tion {kl} leads to a destabilization of its local neighborhood, 
mainly caused by a loss of buttress, but also due to local 
increases in buoyancy and changes in terminus velocity trig- 
gered by calving. In consequence, the stress level in neighbor- 
ing cells {xy} is increased (first term on rhs of (2)). For sim- 



plicity, we assume that the interactions jf, = J(x — k,y — I) 
depend only on the horizontal and vertical distances p = x — k 
and q = y — I between the cells {kl} and {xy}. Further, we 
restrict the model to excitatory (positive) nearest-neighbor 
interactions without self-coupling, i.e. 



J(p, q) = < 



if p = and q = 
if \p\ > 1 or \q\ > 1 
> else . 



(3) 



For the results reported in the next section, we use an asym- 
metric interaction kernel (see Fig. 3C) 



J{p,q)=C{ 



4 if p = and q — 1 

3 if \p\ — 1 and q = 1 

2 if \p\ — 1 and q = 

1 if \p\ < 1 and q = -1 

else . 



(4) 



Here, C is a normalization constant. The asymmetry in the 
vertical direction reflects that cells above the calving cell will 
likely experience a larger stress increment than those below 
due to gravity. To test whether the dynamics of the model 
critically depends on the specific choice of the interaction- 
kernel shape we also consider symmetric kernels 



J{p,q) = C 



1 ifp = and |q| = 1 
I if \p\ = 1 and q = 
else . 



(5) 



Qualitatively, the results for symmetric and asymmetric inter- 
action kernels are the same (not shown here). Note that with 
the symmetric kernel (5), our calving model is (essentially 1 ) 
identical to the sandpile model of Bak et al. (1987, 1988). 

1 ln the model of Bak et al. (1987, 1988), the "stress" z is reset by 
a fixed amount after "calving" , whereas we consider a reset to a 
fixed value 2 = 0. 
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Fig. 3: Sketch of the calving model. A: Time evolution of internal ice stress z in an individual cell. Calving of neighboring 
cells or external perturbations (triangles shown in B) cause jumps in ice stress z. Crossing of critical stress z = 1 (dashed 
horizontal line) leads to calving (triangle-down marker) and reset of stress level to z = 0. C: Glacier front (as seen from the 
sea/fjord; width W , height H) subdivided into WH cells. Calving of cell {kl} (cross) leads to stress increments (gray coded) 
in neighboring cells (depending on relative cell position). 



To study the dependence of the calving dynamics on the 
coupling between cells, we consider the total susceptibility 
w = J(p, q) as a main parameter of the model. It 

characterizes the overall effect of a calving cell on the ice 
stress in its local neighborhood. An increase in the suscepti- 
bility w corresponds to a destabilization of the glacier front. 
Note that w is measured in units of the critical stress z cr jt ; an 
increase in w can therefore also be interpreted as a decrease 
in z cr it- To study the effect of ice susceptibility and/or yield 
stress, it is therefore sufficient to vary w and keep 2 cr it = 1 
constant. Both ice susceptibility and yield stress are deter- 
mined by external factors like temperature, glacier velocity, 
buoyancy, glacier thickness, etc.. An increase in temperature, 
for example, lowers the yield stress (Benn et al., 2007a) and, 
thus, leads to an increase in the susceptibility w. 

Experimental paradigm. Due to the above described 
interactions, calving of individual cells may trigger calving in 
neighboring cells, thereby causing calving avalanches. At the 
beginning of each experiment, the internal ice stress of each 
cell was initialized by a random number drawn from a uniform 
distribution between and 1.1. On average, 10% of the cells 
were therefore above the critical stress z CY [i = 1 and started 
calving immediately. In general, this initial calving stopped 
after some time (see, however, Sec. 3.3). After this warmup 
period, we performed a sequence of perturbation experiments: 
In each trial m — 1, . . . , M, a single cell {kl} was randomly 
chosen and perturbed by a weak delta pulse s^ t (t) = S(t) 
of amplitude J cx t = 0.1 (at the beginning of each trial, time 
was reset to t = 0). The trial was finished when the calving 
activity in response to the perturbation had stopped. We de- 
fine the number of cells calving in a single trial as the event 
size /i. The distance u — v between two subsequent successful 
trials v and it > v, i.e. trials with fi > 0, defines the inter- 
event interval t. Examples of calving activity in individual 
trials are shown in Fig. 6. 



Model parameters. The model parameters are summa- 
rized in Tab. 1. 

Simulation details. The model dynamics was evaluated 
numerically using the neural-network simulator NEST (Gewaltig 
& Diesmann, 2007, see www.nest-initiative.org) which has 
been developed and optimized to simulate large systems of 
pulse-coupled elements. Simulations were performed in dis- 
crete time t = 0, 1, 2, Cell states were updated synchronously, 

i.e. calving activity at time t increments the stress in neigh- 
boring cells at time t + 1. 

2.3. Data analysis 

In the following, we describe the characterization of the marginal 
distributions and auto-correlations of event sizes and inter- 
event intervals. Field and model data were analyzed using 
identical methods. Similarly, we applied identical tools to the 
event sizes /i and the inter-event intervals r. Unless stated 
otherwise, we will therefore not distinguish between fi and r 
in this subsection, and use X as a placeholder. 

Distributions of sizes and intervals. The overall char- 
acteristics of the distribution of data points (i = 1, . . . , n; 
n = sample size) is given by its mean, standard deviation 
SD, minimum and maximum, and the coefficient of varia- 
tion CV = SD/mean (see Tab. 2). For the case of the inter- 
event intervals, the CV provides a measure of the regular- 
ity of the calving process: while CV = corresponds to a 
perfectly regular process with delta-shaped interval distri- 
bution (clock), CV = 1 is characteristic for a process with 
exponential interval distribution, e.g. a Poisson point process 
(Cox , 1962). Histograms of the data on a logarithmic scale 
(relative frequency: number of observations within an inter- 
val [log 10 (X) , log 10 (X) + log 10 (AV)) normalized by n) are 
used for a graphical illustration of the entire distributions. 
As shown in Fig. 5B-G, Fig. 7A,D and Fig. 9 (symbols), the 
empirical distributions obtained this way are broad and re- 
semble power-law or exponential distributions. Note, how- 
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Name 


Description 


Value 


W 


width of glacier front 


{100,200,400} 


H 


height of glacier front 


{25,50, 100} 


-^crit 


critical stress (yield stress) 


1 


W 


calving susceptibility 


{0.5,..., 1.5} 




perturbation amplitude 


0.1 


M 


number of trials 


10000 



Table 1. Model parameters. Curly brackets {. . .} represent parameter ranges. 



ever, that such histograms, obtained by binning of finite data 
sets, are generally biased and therefore not appropriate for a 
quantitative analysis (Clauset et al., 2009). Here, we applied 
maximum-likelihood (ML) methods (see Clauset et al., 2009) 
to quantify to what extent the field and model data X can 
be explained by an exponential or a power-law distribution: 



p cx (X) = iV 0X 
Ppl (X) = N pl 



e- xx 

else 

X ^ < X mm < X < Jfmax 

else 



(6) 
(7) 



The cutoffs were set to the observed minimum and maxi- 
mum, respectively: X min = mhij(JQ), X max = maxj(JSQ). 
The prefactors N p \ and iV ex are normalization constants. The 
exponents A and 7 were obtained by maximizing the log- 



likelihoods I 



cx/pl 



= E,- 



[log(f> cx /pi(X,))J for the two model 
distributions. Here, Ej [. . .] = ^ X/i=l ■ • • denotes the aver- 
age across the ensemble of data points. The quality of the 
ML fits (goodness-of-fit test) was evaluated as described by 
Clauset et al. (2009) using surrogate data and Kolmogorov- 
Smirnov statistics. The resulting p- values indicate how well 
the data can be explained by the model distributions p 0X (X) 
or p p i(X). The log-likelihood ratio R = l p \ — l cx , the differ- 
ence between the maximum log-likelihoods, is used to judge 
which of the two hypotheses, the power-law or the expo- 
nential model, fits the data better. R > indicates that 
the power-law model p p i(X) is superior (and vice versa for 
R < 0). The variance of R, estimated as Ej — R) 2 ] with 
Ri = log(ppi(Xi))— \og(p e x(Xi)) for the best-fit distributions, 
was used to test whether the measured log-likelihood ratio R 
differs significantly from zero (for details, see Clauset et al., 
2009). 

Auto-correlations. To investigate whether calving event 
sizes and intervals are informative about future events, we 
calculated the normalized auto-correlations (see Fig. 11) 



a(i) 



Ej [X(j)X(j + t)] 



(8) 



vnthX j =X j -E j [X,]. 

Software. The data analysis was performed using Python 
(http : / /www . python . org) in combination with NumPy (http : 
//numpy . scipy . org) and SciPy (http : // scipy . org). Results 
were plotted using matplotlib (http: //matplotlib. sourcef orge 
net). 



3. RESULTS 

Calving at the front of tidewater glaciers often occurs as a 
sequence of events, thereby showing the characteristics of 
avalanches: an initial detachment of small ice blocks can cas- 
cade to events of arbitrary size (see Fig. 4). In this article, we 
propose that the underlying dynamics can be understood as a 
result of the mutual interplay between calving and the desta- 
bilization of the local neighborhood of the calving region. By 
means of a simple calving model (see Sec. 2.2), we show that 
this mechanism is sufficient to understand the variability in 
event sizes and inter-event intervals observed in the field data 
(see Sec. 3.1). Fluctuations in external parameters may ad- 
ditionally contribute to the calving variability but are not 
required. This is confirmed by our observation that changes 
in air temperature and tides do not affect the shape of the 
size and interval distributions obtained from the field data 
(see Sec. 3.2). The simple calving model allows us to study 
the effect of glacier parameters on the distributions of event 
sizes and intervals in a controlled manner. An increase in 
the calving susceptibility leads to broader event-size distribu- 
tions. At a critical susceptibility, the model glacier undergoes 
a transition to a regime where a small perturbation leads to 
ongoing calving activity (see Sec. 3.3). Finally, we show that 
the simple calving model is consistent with the field data in 
the sense that the size of future calving events is not corre- 
lated with the size of past events. Predicting event sizes from 
past events is thus difficult, if not impossible (see Sec. 3.4). 

3.1. Variability of event sizes and 
inter-event intervals 

We analyzed data obtained from two glaciers on Svalbard, 
Kronebreen and Sveabreen, during three continuous observa- 
tion periods with more than 7000 calving events in total. The 
longest observation period lasted 12 days with 5868 events 
(Fig. 5A). An overview of the three data sets and the basic 
event-size and interval statistics is provided in Tab. 2. 

The sizes /j, of monitored events are highly variable. They 
extend over 4 orders of magnitudes, from about 10 m 3 up 
to more than 10 5 m 3 (symbols in Fig. 5B-D). The event-size 
coefficient of variation (CV) varies between 6.2 and 8.2; the 
standard deviations are substantially larger than the mean 
values (542-1512 m 3 ). The distributions of event sizes exhibit 
long tails and resemble power laws (7). Maximum-likelihood 
(ML) fitting yields power-law exponents 7^ between 1.7 and 
2 (solid gray lines in Fig. 5B-D). Note, however, that the 
.p-values of the goodness-of-fit test (see Sec. 2.3) are in all 
cases very small, thereby indicating that the power-law model 
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Fig. 4: Typical calving sequence observed on August 16, 2009, 21:46 GMT, at Kronebreen. The detachment of small ice 
blocks (B,C) triggers a large column drop with the entire height of the front collapsing vertically (D-G), followed by large 
column-rotation events with blocks of ice rotating during their fall (H-J). Time distance between consecutive images is 3 s. 
Black arrow and ellipses mark location of individual calving events. 



does not perfectly explain the event sizes. Still, the power-law 
model (7) fits the event sizes fi better than the exponential 
model (6) (dashed curves in Fig. 5B-D); the log- likelihood ra- 
tios R are in all cases significantly greater than zero (R = 2.3 
to 14). 

The inter-event intervals r span more than 2 orders of mag- 
nitude (1 min up to more than 400 min; symbols in Fig. 5E- 
G). The standard deviations exceed the mean durations be- 
tween two events (3-17 min) by a factor of CV = 1.5 to 1.9. 
Hence, the calving process is highly irregular, significantly 
more irregular than a Poisson point process with exponen- 
tial interval distribution (CV = 1; Cox , 1962). The interval 
distributions have a longer tail than predicted by the expo- 
nential model (6) (dashed curves in Fig. 5E-G), but fall off 
more rapidly than the power-law model (7) (solid gray lines in 
Fig. 5E-G). The log-likelihood ratios R w confirm that the 
inter-event intervals r are neither exponentially nor power- 
law distributed. 

The simple calving model can well reproduce the charac- 
teristics of the event-size and interval distributions obtained 
from the field data. An example time series generated by 
the calving model is depicted in Fig. 6A for a susceptibil- 
ity of w = 1.3. Random perturbations of the model glacier 
lead to responses with broadly distributed magnitudes. In 
most cases, there is no response at all (fi = 0). Small events 
are frequently triggered (see example in Fig. 6B,D,F), large 
events are rare (see example in Fig. 6C,E,G). As in the field 
data, the distributions of event sizes fi generated by the calv- 
ing model resemble power-law distributions (Fig. 7A). The 
width of the size distribution increases with the calving sus- 
ceptibility w (Fig. 7B) while the power-law exponent 7 M de- 
creases and approaches 1 for w = 1.3 (solid curve in Fig. 7C). 
Log-likelihood ratios R are always positive and increase with 
w (data not shown), thereby indicating that the power-law 
model (7) fits the data better than the exponential model 
(6). For w — 1.3, the size distribution spans about 6 orders of 
magnitudes (p = 1, ... , 10 6 ). Note that the maximum event 
size can exceed the system size WH (dashed vertical lines in 
Fig. 7A,B) as a single cell can calve several times during one 
event (trial). 



The inter-event intervals r obtained from the calving model 
span about 2 orders of magnitude (Fig. 7D,E). In contrast to 
the event-size distributions, the width of the interval distribu- 
tion is independent of the calving susceptibility w (Fig. 7E). 
ML fitting of the power-law (7) and the exponential model 
(6) yields exponents 7r and A r which are comparable to those 
obtained for the field data (compare Fig. 7F and Fig. 5E-G). 
The log-likelihood ratios R are always close to zero (data not 
shown), again consistent with the field data. 

Note that the calving variability arising from the model is 
not a mere result of the randomness of the external perturba- 
tions (see Sec. 2.2, Experimental paradigm). Restricting the 
repeated external perturbations to one and same cell in the 
center {W/2, H/2} of the front results in slightly narrower 
but qualitatively similar event-size and interval distributions 
(data not shown). 

3.2. Impact of external parameters 

As shown in the previous subsection, the calving model gen- 
erates broad distributions of event sizes and inter-event in- 
tervals, even under perfectly stationary conditions. Fluctu- 
ations in external parameters are therefore not required to 
explain the event-size and interval variability observed in 
the field data. Here, we support this finding by analyzing 
the relation between calving activity and fluctuations in air 
temperature and tides during the observation period. Both 
tides and temperature can, in principle, affect calving activ- 
ity: High tides increase buoyant forces, thereby destabilizing 
the glacier front. An increase in temperature lowers the yield 
stress (Benn et al., 2007a) and therefore leads to an increase 
in the calving susceptibility w. Hence, one may expect that 
temperature and tide-level fluctuations lead to changes in the 
calving-event size and interval statistics, and thereby explain 
the large variability reported in Sec. 3.1 (see Fig. 5). 

Fig. 8 depicts the simultaneous time series for event sizes 
(A), inter-event intervals (B), air temperature (C), change 
in air temperature (D), and tidal amplitude (E) during the 
12-days observation period in 2009 at Kronebreen. Within 
this sampling period, temperatures varied between —0.8 and 
8.8 °C, tide levels between 14 and 178 cm. Significant corre- 
lations between climatic parameters and calving activity are 
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Kronebreen Sveabreen 
2008 2009 2010 



Total number of events 


1041 


5868 


386 


Observation period (days) 


1 


12 


4 


otart date 


26 Aug. 


14 Aug. 


L I J uly 


End date 


1 Sept. 


26 Aug. 


21 July 


Event size fj, (m 3 ) 








Mean 


803 


542 


1512 


Standard deviation SD 


6599 


4014 


9363 


Coefficient of variation CV 


8.2 


7.4 


6.2 


Minimum 


13 


13 


13 


Maximum 


135070 


135070 


135070 


Inter-event interval r (min) 








Mean 


8 


3 


17 


Standard deviation SD 


14 


5 


32 


Coefficient of variation CV 


1.7 


1.5 


1.9 


Minimum 


1 


1 


1 


Maximum 


201 


98 


446 



Table 2. Overview of the three field-data sets with event-size and interval statistics. 



not observed. To test whether fluctuations in air tempera- 
ture and tidal amplitude have an effect on the overall shape 
of the event-size and interval distribution, we grouped the 
data into high/low temperature/tide intervals. The event-size 
and interval histograms obtained for each group are indistin- 
guishable (Fig. 9). Hence, the fluctuations in air temperature 
and tides within the observation period have no effect on the 
shape of the distributions. They cannot explain the observed 
event-size and interval variability. This does not imply that 
changes in temperature and tides do not affect the calving 
statistics in general. Long-term observations are required in 
order to investigate how larger modulations in temperature 
and tide levels (e.g. seasonal fluctuations) affect the shape of 
the event-size and interval distributions (see Sec. 4). 

From the results reported here, we cannot exclude that fluc- 
tuations in other parameters, for example changes in glacier 
velocity or buoyancy, contribute to the size and interval vari- 
ability. Note, however, that both glacier velocity and buoy- 
ancy are affected by calving itself (Benn et al., 2007a). They 
are not purely "external" parameters. In our simple model, 
the positive feedback loop between calving and these param- 
eters is abstracted in the form of the (positive) interaction 
kernel (see Sec. 2.2). Hence, correlations between calving ac- 
tivity and fluctuations in glacier velocity or buoyancy can 
arise naturally from the local calving dynamics at the glacier 
front. 

3.3. Self-sustained calving 

For susceptibilities w < 1.3, small perturbations of the model 
glacier result in calving responses with finite life time (see 
e.g. Fig. 6F,G). Calving of a single cell can trigger an avalanche 



which sooner or later fades away and stops. For w > 1.3, 
the model glacier enters a new regime where the avalanche 
life times seem to diverge (hatched areas in Fig. 7B,C,E,F). 
We illustrate this by randomly initializing the model glacier 
such that a only small fraction (1%) of cells is superthreshold 
(z > 1) at time t — and hence starts calving immediately, 
simulating the model glacier until this initial calving stops 
and measuring the corresponding survival time for a broad 
range of susceptibilities w — 0.8, ...,1.5 (see Fig. 10). For 
w > 1, the survival time quickly increases with w and di- 
verges at about w — 1.3. Beyond this critical susceptibility, 
calving never stopped within the maximum simulation time 
(here 900 time steps). Note that this behavior is robust and 
does not critically depend on the choice of other model pa- 
rameters like the terminus dimensions (W, H) and the shape 
of the interaction kernel (3). In all cases, we observe a critical 
susceptibility close to w = 1.3 (data not shown here). 



3.4. Predictability 

In the field data, the average correlation between the sizes 
fij and fij+i of subsequent events j and j + i hardly differs 
from zero for all i > (Fig. 11A). Hence, predicting future 
event sizes from past events appears hopeless. This behavior 
is reproduced by the simple calving model (Fig. 11B). The 
inter-event intervals r obtained from the field data exhibit 
a moderate long lasting correlation which is not observed in 
the model data (Fig. 11C,D). The reason of this discrepancy 
between the model and the field data is unclear, but we sus- 
pect that it is due to non-stationarities in the field data (see 
Fig. 8B) which are, by construction, absent in the model. 
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Fig. 5: Field data. A: Example time series for 12 observation days at Kronebreen (2009). Each bar represents a calving event of 
size jj, (log scale). B-G: Distributions (log-log scale) of iceberg sizes /i (B-D) and inter-event intervals r (E-G) for Kronebreen 
2008 (B,E), 2009 (C,F) and Sveabreen 2010 (D,G). Field data (symbols), best-fit power-law (solid lines; decay exponents 7^, 
7t) and exponential distributions (dashed curves; decay exponents A^, A T ). R represents corresponding log-likelihood ratio. 



4. DISCUSSION 

We showed that calving-event sizes and inter-event intervals 
obtained from direct continuous observations at the termini 
of two tidewater glaciers are highly variable and broadly dis- 
tributed. We demonstrated that the observed variability can 
be fully explained by the mutual interplay between calving 
and the destabilization of the local neighborhood of the calv- 
ing region. A simple calving model accounting for this inter- 
play reproduces the main characteristics of the event-size and 
interval statistics: i) Event-size distributions resemble power 
laws with long tails spanning several orders of magnitude, ii) 
Interval distributions are broad but fall off more rapidly than 
power laws, iii) Correlations between the sizes of subsequent 
events vanish. We conclude that the observed calving variabil- 
ity is a characteristic feature of calving and is not primarily 
the result of fluctuating external (e.g. climatic) conditions. 
Event sizes of all magnitudes have to be expected, even under 
ideal stationary conditions. 

The calving model predicts that the width of the event-size 
distribution depends on a parameter w representing the calv- 
ing susceptibility of the glacier terminus: Roughly, w mea- 
sures how prone the glacier is to calve in response to calv- 
ing. It describes to what extent calving increases the internal 



ice stress in the neighborhood of the calving region. Alterna- 
tively, w may represent the inverse of the yield stress, i.e. the 
critical stress at which ice breaks. The susceptibility w is 
determined by the properties of the ice and by factors like 
temperature, glacier velocity, buoyancy or glacier thickness. 
In our model, the width of the event-size distribution in- 
creases monotonously with w. Therefore, the shape of the size 
distribution, as-for example-characterized by the power-law 
exponent 7^, may be informative about the stability of the 
glacier terminus. In contrast to the event-size distributions, 
the shape of the inter-event interval distribution is insensi- 
tive to the susceptibility w. At a critical susceptibility w CI n, 
the model predicts an abrupt transition of the glacier to a 
new regime characterized by ongoing, self-sustained calving 
activity. Observations of rapid glacier retreats (Pfeffer, 2007; 
Briner et al., 2009; Motyka et al., 2011) may be explained 
by these supercritical dynamics (see also Amundson & Truf- 
fer, 2010). In the model, the power-law exponent 7^ is close 
to 1 as the susceptibility approaches the critical value w CT n 
(see solid curve Fig. 7C). We found that this observation does 
not critically depend on the choice of the model parameters 
(glacier dimensions, shape of the interaction kernel, pertur- 
bation paradigm). This suggests that the power-law expo- 
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Fig. 6: Size and duration variability of calving events in response to random external perturbations (model results). A: Glacier 
responses (event size n; log scale) to 1000 consecutive random perturbations. Black square (trial 190) and circle (trial 884) 
mark events shown in B,D,F and C,E,G, respectively. B-G: Moderate-size (B,D,F; /i = 15) and large calving event (C,E,G: 
H — 27903) in response to punctual random perturbations (perturbation sites indicated by black stars) in trials 190 and 884, 
respectively. B,C: Distribution of stress z X y (gray coded) across the model-glacier front after calving response. D,E: Spatial 
spread of single-trial calving activity across the model glacier front. Gray dots mark cells which have calved. F,G: Spatio- 
temporal spread of the events shown in D,E (id of cell at position (x,y) = Hx + y; discrete time). Note different scales in F 
and G. Glacier width W = 200, glacier height H = 50, calving susceptibility w = 1.3. 



nent 7^ may serve as an indicator of a glacier's proximity to 
the transition point where it starts retreating rapidly. Bassis 
(2010) proposed a similar stability criterion which depends 
on various geometric and dynamical near-terminus parame- 
ters. In our case, the diagnostics is exclusively based on the 
distribution of event sizes. 

The calving model in this study is inspired by previous 
work on the emergence of power-law distributions. Power- 
law shaped magnitude distributions are abundant in nature. 
They are found, for example, for earthquake magnitudes (Bak 
et al., 2002; Hainzl, 2003), the luminosity of stars (Bak, 1996), 
avalanche sizes in sandpiles (Bak, 1996), landslide areas (Guzzetti 
et al., 2002), subglacial water pressure pulses (Kavanaugh, 
2009), dislocation avalanches in ice (Richeton et al., 2005), 
and sea ice fracturing (Rampal et al., 2008). Bak et al. (1987) 
demonstrated that power-law distributions can arise natu- 
rally in spatially extended dynamical systems which have 
evolved into a "self-organized critical state" consisting of min- 



imally stable clusters of all length scales. Perturbations can 
propagate through the system and evoke responses character- 
ized by the absence of spatial and temporal scales (avalanches). 
The calving model used in the present study is qualitatively 
similar to the sandpile model of Bak et al. (1987). It is there- 
fore not surprising that it predicts power-law shaped event- 
size distributions. The scientific value of this part of our work 
exists in mapping the original model by Bak et al. (1987) to 
the dynamics of glacier calving and in relating the model pa- 
rameters to physical measures. Moreover, the phenomenon of 
self-sustained ongoing activity has, to our knowledge, never 
been reported so far in this context. 

Several aspects of this work may be subject to criticism 
and improvements, in particular the data acquisition and the 
construction of the calving model: Direct visual observation 
by humans is so far one of the most reliable methods to mon- 
itor individual calving events in continuous time (van der 
Veen, 1997). The results are however hard to reproduce ex- 
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Fig. 7: Size (top row) and inter-event interval distributions (bottom row) in the calving model. A,D: Distributions (log-log 
scale) of iceberg sizes fi (A) and inter-event intervals r (D) for calving susceptibility w — 1.25. Simulation results (symbols), 
best-fit power-law (solid lines) and exponential distributions (dashed lines). B,E: Dependence of size and interval distributions 
(log-log scale) on calving susceptibility w. White curves mark maximum sizes and inter-event intervals. C,F: Dependence of 
decay exponents of best-fit power- law (solid lines; 7^, 7 T ) and exponential distributions (dotted lines; X^, X T ). Glacier width 
W — 400, glacier height H = 100. Vertical dashed lines in A and B indicate system size WH — 40000. Dotted horizontal and 
vertical lines in B,E and C,F, respectively, mark susceptibility w — 1.25 used in A and D. Hatched areas in B,C,E,F correspond 
to regions with ongoing, self-sustained calving (see Sec. 3.3). 



actly: Due to the sporadic unpredictable nature of calving 
it is difficult to capture each event. Different observers may 
be in different attentive states. Each new observer needs to 
adjust her- or himself to a common scale of perceived sizes 
yj. To minimize the variability in perceived sizes we arranged 
test observations where all observers were simultaneously con- 
fronted with size estimates (see Sec. 2.1). To confirm our find- 
ings on the relation between the shape of the event-size dis- 
tribution and the state of the glacier, more data are needed 
from glaciers in different environments (freshwater, floating 
tongue, iceshelves) and different dynamical states (advanc- 
ing, stable, retreating) . In order to test our stability criterion 
one would need to continuously monitor the same glacier year 
after year for a few days. Long-term continuous observations 
or observations in cold months are difficult. Automatic mon- 
itoring would therefore be highly desirable. Unfortunately, 
all available automatic monitoring methods have limitations. 
Terrestrial photogrammetry, for example, is limited by the 
iceberg size, visibility and illumination of the glacier. Ice- 
berg size and type also limit the use of ground-based radar. 
Chapuis et al. (2010) showed that radar could detect only 
events larger than 150 m 3 . Remote sensing (optical and radar 
imagery) has the same limitations as terrestrial photogram- 
metry. In addition, its low temporal resolution does not al- 
low to register individual calving events. Seismic monitoring 
(O'Neel et al., 2010) is a very promising technique but can 



detect only the largest events (Kohler et al., 2012). The range 
of event sizes accessible by seismic methods may however be 
sufficient to estimate power-law exponents. A major prob- 
lem in studying the statistics of individual calving events is 
to define what a single event actually is. In the framework 
of our model, an event is defined as the total response to 
a single perturbation. In our numerical experiments, a new 
perturbation is not applied before the response to the previ- 
ous perturbation has stopped. In nature, however, the glacier 
front may be constantly perturbed, e.g. by the movement of 
the glacier. Several events may be triggered simultaneously 
in neighboring regions of the front and, hence, overlap both 
spatially and temporally. The separation of individual events 
in the field data can therefore be difficult. In consequence, 
short intervals and small events may be underestimated. 

The calving model used in this study is highly minimalistic. 
It implements the mechanism which we think is essential for 
an understanding of the calving variability observed in the 
field data (i.e. a positive feedback loop between calving and 
destabilization of the glacier front), but neglects a variety of 
factors: First, the model is two-dimensional. Stress can only 
propagate tangentially and not in the third dimension per- 
pendicular to the front. Extending the model to three dimen- 
sions is not straightforward, as it is unclear how calving at 
the front affects regions deeper within the ice body. Second, 
the interaction kernel is restricted to nearest-neighbor inter- 
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Fig. 8: Relation between calving and external parameters (field data, Kronebreen, 2009). A: Event sizes fi (log scale). Individual 
events (dots) and 4-hour average (black curve). B: Inter-event intervals r (log scale). Individual events (dots) and 8-hour average 
(black curve). C: Air temperature (Ny-Alesund, Norwegian Meteorological Institute). D: Change in air temperature within 6 
hours. E: Tidal amplitude (Norwegian Mapping Authorities). 



actions. Whether and how calving affects regions more distant 
from the calving region is unclear. A direct measurement of 
the interaction kernel in the field is difficult as it would require 
monitoring of changes in ice stress in response to individual 
calving events. Conclusions on the shape of the interaction 
kernel could be drawn indirectly from the spatial structure of 
calving avalanches (e.g. obtained by terrestrial photogramme- 
try). For illustration, consider Fig. 6E: The triangular shape 
of this large event is a direct consequence of the asymmetry of 
the interaction kernel (4) (see Fig. 3C). For symmetric kernels 
(5), we observed that the spatial structure of calving events is, 
on average, symmetric (data not shown here). Note, however, 
that for the main findings of our study, the exact shape of the 
interaction kernel is not critical. Third, the relations between 
the susceptibility w (the area of the interaction kernel, see 
Sec. 2.2) and external parameters like air temperature, tides, 
glacier velocity or buoyancy are unclear. These relations could 
be established empirically from a larger data set obtained 
from continuous long-term observations of different glaciers 
during different dynamical states. For each identified set of 
external parameters, our model could be fitted to the distribu- 
tion of monitored event sizes, thereby providing an estimate 
of w. Based on these relationships, the event-size distributions 
for a new "test" data set (i.e. for data not used for the fitting 
procedure) could be predicted by the model and compared 
to the field-data distributions. Fourth, stress increments in 



response to calving or external perturbations are instanta- 
neous in the model. In reality, these interactions are likely to 
be smoother and delayed. Detailed information about this is 
however limited. Fifth, in the absence of calving and exter- 
nal perturbations, internal ice stress in a given model cell is 
constant. In real ice, stress in a certain volume element may 
slowly dissipate and decay to zero. We assume that the time 
constant of this decay is much larger than the time between 
two (external or internal) perturbations and therefore can be 
approximated as being infinite. Sixth, external perturbations 
are modeled as punctual (delta-shaped) events in time and 
space to characterize the system's response in a well defined 
manner (i.e. to evoke responses without spatio-temporal over- 
lap; see above). In reality, external perturbations are spatially 
and temporally distributed (e.g. glacier movement, glacier ve- 
locity gradients, buoyancy). Finally, the model neglects sub- 
marine calving which, in reality, amounts for about 13% of 
the total ice loss at the front. One way to account for the 
submarine dynamics would be to assign different calving sus- 
ceptibilities to the subaerial and the submarine parts. Our 
model describes calving at the level of individual events and 
focuses on a single aspect of calving: the interplay between 
calving and front destabilization. Factors like air tempera- 
ture, water depth, height-above-buoyancy, surface melt, ice 
thickness, thickness gradient, strain rate, mass-balance rate, 
backward melting of the terminus, etc., are included indi- 
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Fig. 9: Distributions of iceberg sizes fi (A) and inter-event intervals r (B) for four different combinations of air-temperature 
and tide ranges (Kronebreen, 2009; log-log scaling): high temperature - high tide (HH), high temperature - low tide (HL), 
low temperature - high tide (LH), low temperature - low tide (LL). Low/High air-temperature interval: — 0.8-4.2/4.2-8.8 °C, 
Low/High tide-level interval: 14-92/92-178 cm. 



rectly or lumped together and treated as constant parame- 
ters. Previous macroscopic models (Brown et al., 1982; van 
der Veen, 1996; Vieli et al., 2001, 2002; Benn et al., 2007b; 
Amundson & Truffer, 2010; Nick et al., 2010; Otero et al., 
2010; Oerlemans et al., 2011) relate these factors to the over- 
all calving rate. A description of size and timing of individual 
calving events is beyond the scope of most of these mod- 
els. Our model could be connected to these models by link- 
ing the overall calving rate to the calving susceptibility or 
the characteristics of external perturbations. The resulting 
multi-scale model would have the potential to describe the 
statistics of event sizes and inter-event intervals under more 
realistic non-stationary conditions. The rationale behind the 
simplistic approach presented in our study is the opposite: 
to construct a minimal model which reproduces the main 
characteristics of the event-size and interval statistics. Dis- 
connecting the glacier front from the external world allows us 
to study calving dynamics under ideal stationary conditions 
and to demonstrate that the observed calving variability is 
not primarily a result of fluctuations in external conditions 
but can emerge from simple principles governing the dynam- 
ics directly at the glacier front. 
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